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Control of Superconducting Correlations in High-Tc Cuprates 
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A strategy to enhance d-wave superconducting correlations is proposed based on our numer- 
ical study for correlated electron models for high-Tc cuprates. We observe that the pairing is 
enhanced when the single-electron level around (ir, 0) is close to the Fermi level Ef, while the 
d-wave pairing interaction itself contains elements to disfavor the pairing due to shift of the 
(n, 0)-level. Angle-resolved photoemission results in the cuprates are consistently explained in 
the presence of the d-wave pairing interaction. Our proposal is the tuning of the (tt, 0)-level 
under the many-body effects to Ef by optimal design of band structure. 

KEYWORDS: high-Tc superconductivity, pairing correlations, single-electron level, angle-resolved photoemission 



We report a new aspect of control of superconducting 
correlations in high-Tc cuprates. Our new observation 
in the light of numerical studies is about the critical in- 
terplay between single-electron level at wave vector (tt,0) 
and (0,7r) and electron-pair hopping processes, which are 
usually overlooked. One potential implication is that the 
superconducting transition temperature T c may be con- 
trolled by tuning the single-electron level and the pair- 
hopping coupling constant in a self-consistent fashion. 
This could provide a guideline for the quest of new ma- 
terials with higher T c . 

Since Jjhe first report of high-Tc superconductivity 
in 1986^ various materials with the common Cu02 
layer structure have been found to show superconduc- 
tivity with T c up to 133 [K-.in mercury compounds' 
and 164 [K] under pressured One noticeable char- 
acter common among the high-Tc family is the pres- 
ence of a flat band around the Of.0) and equivalent 
(0,7r) points in the Brillouin zonccP The anomalous 
flat dispersion determined by the angular-resolved pho- 
toemission (ARPES) experiments may be fitted by a 
fc 4 -type formO-* and numerical calculations on the two- 
dimensional (2D) Hubbard and t-J models indicate that 
it is due to strong, glectron-electron correlations rather 
than band effects .§0BB> ARPES data in fact show the 
d-wave gap starts growing feora .the. flak-dispersion re- 
gion in underdoped samplesJilJ'OBEj'O suggesting a 
key role of this region. However, in the insulating or 
very lightly doped samples, the broad quasiparticle peak 
at (ir,0) appear^ substantially lower than the level at 
(7r/2, 7r/2)Ji3'HE3E-!P which implies rapid change in the 
shape of the Fermi surface with increasing doping. Ex- 
perimentally, the rigid band picture thus fails and it is 
yet clear how the electron correlation and pairing inter- 
action effects modify the single-particle dispersion. 

In this letter, we first show how the dynamical su- 
perconducting correlations are enhanced by controlling 
the single-electron level at the wave vector k = (n,0), 
and we will show later that self-consistent determina- 
tion of the flat-dispersion level around (tt, 0) under strong 



correlation effects opens a way to contmL-Lhe. enhance- 
ment. We employ the 2D t-J modei§BE3) for our 
numerical investigation, and introduce the second and 
third neighbor hoppings in order to control this en- 
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where the hopping integrals are ty = t for the nearest- 
neighbor pairs, t\ and ti for the second and third neigh- 
bor pairs, respectively, and otherwise 0. The values 
of t and J are determined from various experiments, 
and a good estimate is t=0.3[eV] and J=0.12[eV] in all 
high-Tc compounds. Hereafter, we will use these val- 
ues unless mention explicitly, and measure energy in 
units of t. The single-electron kinetic energy is then 
given by £n(k)=— 2i(cos k x + cos k y ) — 4tj cos k x cos k y — 
2i2(cos2fca; + cos2fcj,), and the tight-binding Fermi sur- 
face changes its shape at half filling from the 7r/4-rotated 
square due to the added t\- and i2-terms. Fermi surface 
in various high-Tc corooounds has been determined by 
ARPES experiments^^ and based on comparison with 
the hand calculation: for example, for YBa2Cu307, see 
Ref.Eif A couple of previous studies argue that these 
distant -hrxppings alone account for the change in Fermi 
surfaced EJf Another explanation is provided for the 
deep (xpO)-level, based on the spin-charge separation sce- 
narioBJ' We take an alternative interpretation on this 
difference to account for superconducting correlations. 

In order to investigate the enhancement of supercon- 
ductivity, we have calculated its dynamic correlation 
function, Pd(u), instead of the usual equal-time corre- 
lations. It is essential to separate the coherent part to 
measure the enhancement, since only the low-energy dy- 
namics of Cooper pairs is relevant to the superconduct- 
ing transition. Typical results are shown in Figure 1 
for various t\ and ti for creating a spin-singlet hole pair 
on neighboring sites, A = J2j a /rf( a ) c j'T c j+al with the 
d x 2_ y 2 form factor fd(a) = +1 (—1) when a = ±x (±y). 
This is calculated at zero temperature from the ground- 
state wave function obtained by exact diagonalization for 
the 4x4-site cluster. The cluster is chosen such that we 
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Fig. 1. Dynamic correlations of electron-pair annihilation (solid) 
and creation (dotted) with the d x i_ y i symmetry in the 2D t-J 
model. J/t=0A. 



can compare the electron level between the two charac- 
teristic k-points on the Fermi surface at half filling, (ir,0) 
and (7r/2,7r/2). Note that in the noninteracting system, 
the difference in kinetic energy, e (ir, 0) — e (n/2, ir/2), 
is given by uq = 4(ii — 2£ 2 )- 

An important result is that upon doping a single 
hole pair into the half-filled insulator the coherent peak 
Pd(oj — 0) grows monotonically with increasing uq. This 
is a clear evidence of the evolution of coherence of hole- 
pair motion. On the other hand, when uo is negative and 
large, the coherent peak loses its weight, which means the 
motion of the hole pair is substantially damped. We can 
easily understand this behavior by considering the single- 
electron level. For larger u , s (k) is higher around (n, 0) 
and the energy cost of creating a hole pair is smaller 
there. The c4.2_ y 2-wave Cooper pairs under consider- 
ation consist mainly of those holes, and consequently 
the superconducting correlations are enhanced. The Mo- 
dependence is no longer monotonic when a hole pair is 
added to the systems already doped, as shown in Figure 
1 for the doping 5=0.125, but a similar interpretation 
holds. For un > 0, electrons have been already removed 
from the (tt, 0) and (0,7r) points at finite dopings, and 
an additional hole pair can hardly reside around these 
k-points. On the other hand, when uq < 0, electrons 
are removed first from the (±7r/2, ±7r/2) and then from 
(tt, 0) and (0, 7r) upon hole doping. Therefore, when uq is 
negative, the electron occupancy at a finite doping may 
be such that the (it, 0)-level is still below and close to the 
Fermi energy. When this is the case, the superconduct- 
ing correlations will be enhanced. 

Since the above discussion on the single-electron level 
is within the free electron picture, we have to check it di- 
rectly through the single-electron Green's function. For 
the same set of parameters as in Fig. 1, the spectral func- 
tion of the Green's function, A-^(uj), is calculated and the 
results are plotted in Fig. 2. Fermi energy is here set to 
be zero, Ep=0. Sharp peaks near Ep are of quasipar- 
ticle or quasihole, and a broad incoherent background 
is spread at higher \lu\. Here, we define the difference 
in the quasiparticle level u = s(tt, 0) — e(7r/2, ir/2), from 
the peak positions in (to) . The observed w-dependence 
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Fig. 2. Single-electron spectral function in the 2D t-J model in- 
cluding t± and t2- J/t = 0.4. Calculations are performed for all 
the six nonequivalent k-points, and the results for typical values 
of t± and t2 are shown here for k=(7r, 0) and (tt/2, n/2). Parti- 
cle and hole parts correspond to uj > and u> < 0, respectively. 
Fermi energy at <5 = is set that the highest hole peak is at 
w = 0. 



qualitatively support the picture explained before. The 
quasiparticle level at (tt,0) shifts downwards with de- 
creasing wo, despite the reduction of the size of the shift. 
The general trend holds that the (7r,0)-level is close to 
Ep when the pairing correlations are enhanced. Similar 
behavior is also observed at other J /Vs. The reduction 
of the shift is due to electron-electron correlations, since 
single-electron hopping processes lose dispersion and co- 
herence near half filling. We note that the t-J model with 
ti=*2=0 gives u ~ 0, and this degeneracy alscLxemains 
in much larger systems of the Hubbard model E3 
As is mentioned above, the ARPES results^ 
suggest that the (w, 0)-level is substantially lower than 
the (tt/2, 7r/2)-level in the insulating or very lightly 
doped samples. From the above argument, if we employ 
the t-J model, the deep (w, 0)-level requipes_that the dis- 
tant hopping satisfies uq < as in RefsEj'H However, 
if this is the case, the pairing correlations are not op- 
timally enhanced at realistic dopings around 6 ~ 0.15. 
Experimentally, the (jr, 0)-level quickly rises up to the 
(tt/2, 7r/2)-level with further doping. This behavior is 
difficult to explain in the rigid band picture, evidence of 
correlation effects presumably including strong pairing 
interaction. 

To incorporate such pairing effects, we now examine 
microscopic processes which account for generic feature 
of d-wajie pairing by employing the square of local kinetic 
energy .Ef A part oi^linilaj: interaction has been studied 
by several authors ,E3E3Ej'E$ but many features remain 
unexplored. The Hamiltonian, say the t-J-Wd model, 
reads 

H tJWd =H t j-W d Y,^fd{a)c) a c j+a „+tt.c]\ (1) 

where a labels four nearest neighbors of each site, and the 
form factor fd(a) is d x 2 _ y 2-w&ve like as defined before. 
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Fig. 3. Two typical third-order terms in the strong coupling ex- 
pansion of the d-p model leading to the W<j term. Numbers 1-3 
denote the order of the perturbation processes, and the process 
2 uses the direct O-O hopping. Arrows show hole spins, and the 
Zhang- Rice singlets are depicted by dotted ellipses, (a): Three- 
site term with 90° configuration, (6): 180° configuration. Note 
that the number of equivalent intermediate configurations differs 
between (a) and (6). 



The VKj-jterm consists of two-site terms and three-site 
terms. ej'EZP The former renormalizes the superexchange 
coupling constant as J e s=J + 8Wd, and its effects could 
be already contained in those of the J-term. On the con- 
trary, the three-site terms have new effects. Since it may 
be rewritten into the form of electron-pair hoppings ,EZf 
it is natural to expect further enhancement of super- 
conductivity due to these terms. This W^-term may be 
derived as the low-energy effective Hamiltonian starting 
from the_d-p modeljL!! as Zhang- Rice construction of the 
J-termJ13) but has customarily been omitted without ex- 
amining its relevance. One may also obtain the three-site 
term starting from the single-band -Hubbard model in 
the limit of strong on-site repulsion pS but then with a 
constant form factor, /(£) — » Lpandits effects have been 
studied by several authors.E^L^ElJ'EP Generally, starting 
from the d-p model, the W-term with the d x 2_ y 2-wave 
form factor, Wd, and that with the constant form factor, 
W s , are both obtained. Figure 3 sketches two of typi- 
cal fundamental processes contributing to the three-site 
terms. The difference of the two contributes to the W-- 
term, while the sum leads to the W^-term. Although 
usually neglected, these processes are in fact of lower 
order than those for superexchange, meaning a consider- 
able size of the coupling constants, Wd and W s . 

In the following we will demonstrate based on numer- 
ical calculations an important role of the W^-term for 
the enhancement of superconducting correlations. We 
also show that such pairing interaction simultaneously 
induces renormalization of single-particle dispersion. In 
contrast to the W^-term, the W s -term does not affect 
the single-electron level at k's along the Fermi surface at 
half filling. c3 Therefore, the deep (7r,0)-level observed in 
experiments again require uo < if Wd — and W s ^ 0, 
and this implies that the enhancement of superconduc- 
tivity remains not prominent. On the other hand, the 
Wd-term drives pairing instability, and simultaneously 
shifts the single-particle (7r,0)-level downwards relative 
to (7r/2,7r/2) as we show below. It turns out that these 
two phenomena are compatible only when the W^-term 
is playing the driving mechanism of pairing. Based on 
these observations, we examine how to optimize super- 
conductivity with the VFd-term. 

Figure 4 shows the dynamic superconducting correla- 
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Fig. 4. Dynamic correlations of electron-pair annihilation in the 
2D t-J-Wd model at 6 = 0. ti = t 2 = and J/t = 0.3. 

tions of the t-J-Wd model (|l]), and here t\ = ti = 0. 
The growth of the coherent peak Pd{uj) upon switching 
Wd is noticeable, and at the same time the weight of 
the incoherent background at higher energies is substan- 
tially suppressed. Thus the effects of the M^-term on the 
enhancement of superconducting correlations are clear. 

The W^-term modifies single-electron dynamics as 
well. Electrons hop using this term to second and third 
neighbor sites with and without spin flip accompanied, 
where the single-electron part of the W^-term renormal- 
izes these hoppings as ti + 2Wd and t^ — Wd, respectively. 
Quasiparticle energies e(k) are affected accordingly. Of 
course, electrons are strongly correlated near half filling, 
and we need numerical calculations to determine the real 
quasiparticle energy. We have calculated the Green's 
function for the t-J-W^ model now including t\ and ti, 
and determined the hole occupation and the quasiparti- 
cle/quasihole energy at each k. Some results are shown 
in Fig. 5. As was shown before, the (5-dependence of the 
single-electron level is basically consistent with the rigid 
band picture, when Wd — 0. On the other hand, for finite 
Wd, the single-electron level is strongly renormalized to 
gain a large energy from the W^-term rather than single- 
electron hopping term. Its important consequence is the 
pinning of the (7r,0)-level near Ep during hole doping, 
and the spectral function is found to have a sharp low- 
energy peak in both the particle and hole parts. This 
double peak structure agrees with what is expected in 
superconducting states, and it is more prominent than 
that in Fig. 2. 

The single-electron part of the W^-term, to some ex- 
tent, suppresses superconductivity near half filling, since 
it pushes down the (7r,0)-level away from Ep. On more 
general grounds, the d-wave pairing interaction neces- 
sarily induces a downward shift of (7r,0)-level. The ex- 
perimentally observed {ir,0) level lower than (ir/2,ir/2) 
in the highly underdoped samples may at least partially 
result from such mechanism. We have calculated pair- 
ing correlations for the t-J-W^ model including t\ and 
ti hoppings and examined the relation with the (7r,0)- 
level. Figure 6 shows its coherent part Pd(u = 0) at 
S = as a function of the level difference between (tt,0) 
and (7r/2,7r/2). As a realistic choice, we here set J = 
and Wd — 0.05i, equivalent to J c s = OAt, corresponding 
to Fig. 1. The level difference u < at half filling and 
the enhancement of the pairing correlations with further 
doping now becomes more compatible. In addition, as in 
the t-J model, superconductivity is more enhanced when 
the (7r,0)-level is tuned close to the Fermi energy, and 
this indeed offers potential optimization beyond present 
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Fig. 6. Enhancement of superconductivity by level tuning in the 
2D t-Wd model. The coherent part of electron-pair annihilation 
at 5 = is plotted as a function of the level difference u between 
(7r,0) and (7r/2, 7r/2). Level difference for the free electron case 
uq is also shown for comparison. 



Fig. 5. Single-electron spectral function in the 2D t-Wj model 
including t\ and t2- J/t = and Wd/t = 0.05. Results for k= 
(tt,0) and (7r/2,7r/2) are shown here for typical values of t\ and 
t 2 - 



experimental realizations. 

We briefly note an important role of the dynamic ex- 
ponent z of the Mott transition. It sets the coherence 
temperature as T* ~ <5 z / 2 j33>£l) an essential factor deter- 
mining T c . Scaling of superfluid density is also subject to 
z, since its upper bound is given by Drude weightt9 and 
this scales as D ~ 8 Z / 2 in 2D as 5 -> 0B> Therefore, for 
larger z, the superfluid density vanishes faster with ap- 
proaching half filling. When the W s -term is switched 
on, we have found the change in z from the unusual 
value z=A to 2, Ex indicating an incoherent-to-coherent 
transition. We speculate that the enhancement of super- 
conductivity we observed near half filling is also associ- 
ated with this transition of z caused by the presence of 
the Wd-term together with the flat band around (71", 0) 
pinned near Ep. We will discuss this point in more detail 
in a future publication. 

In conclusion, the flat dispersion near (it, 0) is deter- 
mined self-consistently from strong correlations and pair- 
ing effects, and the d-wave superconducting correlation 
is enhanced when the resultant band is pinned near the 
Fermi level. Because of its large dependence, it opens 
a possibility of further optimization of pairing correla- 
tion, and hence a chance to achieve a higher transition 
temperature by careful tuning. In the light of this level 
tuning, the high-Tc cuprates may not be in the optimized 
condition, because of the apparently lower level at (71", 0) 
than (7r/2, 7r/2). This is presumably due to the same 
origin with the pairing force, and there is a good chance 
for further tuning. Our proposal is that the bare single- 
particle transfer should be chosen to keep the renormal- 
ized flat dispersion near the Fermi level by compensating 
the downward shift of the (7r,0)-level induced by the d- 
wave pairing interaction itself. This offers a typical way 
to achieve the enhancement of superconductivity by a 
careful tuning of the lattice parameters. For the quest 
of new materials, the bare single-particle level suggested 
from the band structure calculation would be helpful if 
combined with our procedure. 



The authors thank Fakher Assaad for valuable dis- 
cussions. This work is supported by a Grant-in-Aid for 
"Research for the Future" Program from Japan Society 
for the Promotion of Science under the project JSPS- 
RFTF97P01103. 



[7; 

[8 
[9 
[10 

[11 
[12 
[13 
[14 



[15. 
[16 
[17[ 

[18 
[19 
[20 

[21 
[22 

[23 
[24 
[25 
[26 
[27 
[28 
[29 

[30 

[31 

[32 
[33 



J. G. Bednorz and K. A. Miiller: Z. Phys. B 64 (1986) 189. 

A. Schilling et al.: Nature 362 (1993) 56. 
L. Gao et al: phys. Rev. B 50 (1994) 4260. 

K. Gofron et al.: Phys. Rev. Lett. 73 (1994) 3302; Z. X. Shen 
and D. S. Dessau: Phys. Rep. 253 (1995) 1. 
T. Takahashi: private communication. 

M. Imada, et al.: Physics and Chemistry of Transition Metal 
Oxides, ed. by H. Fukuyama and N. Nagaosa (Springer Verlag, 
Berlin, 1999), p. 120. 

F. Assaad and M. Imada: Europhys. J. B in press. 
E. Dagotto et al.: Phys. Rev. Lett. 73 (1994) 728. 
N. Bulut et al.: Phys. Rev. B 50 (1994) 7215. 
D. S. Marshall ct al.: Phys. Rev. Lett. 76 (1996) 4841. 
H. Ding et al.: Nature 382 (1996) 51. 

D. M. King et al.: J. Phys. Chem. Solids 56 (1995) 1865. 
M. R. Norman et al.: Nature 392 (1998) 157. 
J. C. Campuzano et al.: Physics and Chemistry of Transition 
Metal Oxides, ed. by H. Fukuyama and N. Nagaosa (Springer 
Verlag, Berlin, 1999), p. 152. 

B. O. Wells et al.: Phys. Rev. Lett. 74 (1995) 964. 

C. Kim et al.: Phys. Rev. Lett. 80 (1998) 4245 



A. Ino et al.: cond-mat/9809311 (1998); cond-mat/9902048 
(1999). 

P. W. Anderson: Science 235 (1987) 1196. 
F. C. Zhang and T. M. Rice: Phys. Rev. B 37 (1988) 3759. 
See for a review of numerical results, E. Dagotto: Rev. Mod. 
Phys. 66 (1994) 763. 

M. Imada et al.: Rev. Mod. Phys. 70 (1998) 1039. 
R. Lin ct al.: Phys. Rev. B 46 (1992) 11056; W. E. Pickett 
et al.: Phys. Rev. B 42 (1990) 8764. 
A. Nazarenko et al.: Phys. Rev. B 51 (1995) 8676. 
B. Laughlin: Phys. Rev. Lett. 79 (1997) 1726. 
Assaad and M. Imada: unpublished. 
Assaad and M. Imada: Phys. Rev. B 58 (1998) 1845. 
Assaad et al.: Phys. Rev. B 56 (1997) 15001. 
E. Hirsch: Phys. Rev. Lett. 54 (1985) 1317. 
Ammon et al.: Phys. Rev. B 52 (1995) 629, and references 
therein. 

H. Tsunetsugu and M. Imada: J. Phys. Soc. Jpn. 67 (1998) 
1864. 

H. Tsunetsugu and M. Imada: unpublished. 

M. Imada: J. Phys. Soc. Jpn. 64 (1995) 2954. 

D. J. Scalapino et al.: Phys. Rev. B 47 (1993) 7995. 



